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A truncation scheme for the Dyson-Schwinger equations of quantum chromodynamics in Coulomb 
gauge within the first order formalism is presented. The truncation is based on an Ansatz for the 
Coulomb kernel occurring in the action. Results at leading loop order and in the infrared are 
discussed for both the Yang-Mills and quark sectors. It is found that the resulting equations for the 
static gluon and quark propagators agree with those derived in a quasi-particle approximation to the 
canonical Hamiltonian approach. Moreover, a connection to the heavy quark limit is established. 
The equations are analyzed numerically and it is seen that in both the gluonic and quark sectors, a 
nontrivial dynamical infrared mass scale emerges. 
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O . I- INTRODUCTION 

Confinement and the dynamical breaking of chiral symmetry are two highly nontrivial, nonperturbative aspects of 
^vq ' the hadron spectrum. Ideally, we would like to understand both from first principles calculations of the underlying 
gauge theory which is quantum chromodynamics (QCD). Because both effects are manifested in the infrared regime 
(and where singularities may occur), nonperturbative continuum functional methods are one framework within which 
to study the problem. In such studies, one is invariably forced to choose a gauge and our choice is Coulomb gauge, 
r^' One initial reason for this choice is that there exists an appealing picture for confinement: the Gribov-Zwanziger 
scenario [1-3]. 

^ ■ Coulomb gauge studies of QCD have become an area of increasing importance in recent years, mainly due to 
progress within the canonical Hamiltonian approach (see, for example, Refs. [4-11] and references therein), although 
such studies have been around for quite some time (e.g., Refs. [12, 13]). The essential idea is that given the Hamilton 
operator [14] and an Ansatz for the ground state vacuum wavefunctional (supplemented with Gauss' law), the vari- 
ational principle can then be used to generate equations for the various Green's functions of the theory. For a given 
Ansatz, the resulting equations are exact. 
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A second approach to nonperturbative QCD in Coulomb gauge (and that considered here) is to study the Dyson- 
\^ Schwinger equations [15]. In contrast to the canonical approach, Dyson-Schwinger studies in Coulomb gauge are far 
• ■ less developed mainly due to their inherent technical challenges (although attempts have been made [16, 17]). This 
I . being said, various facets of the formalism have been understood: their one-loop perturbative behavior [18-20], the 
' Slavnov- Taylor identities [21], the emergence of a nonperturbative constraint on the total charge [22] and the case 
. of heavy quarks [23-25]. The standard Dyson-Schwinger approach is distinct from the canonical approach in two 
IL^ ' primary respects: it includes the energy dependence of the Green's functions from the outset and must therefore deal 
directly with the noncovariance of Coulomb gauge; also, truncations are typically made to the individual terms in the 
(full) equations. 

] Clearly, it is desirable to be able to connect the canonical Hamiltonian and Dyson-Schwinger formalisms. Working 
in concert, the two different formulations can mutually reinforce each other - what is difficult in one may be intuitive 
in the other and vice versa. An initial aim of this paper is thus to show how the gap equations for the static gluon 
and quark propagators, obtained originally within the canonical formalism [6, 12], can be derived from a leading 
order truncation of the Dyson-Schwinger equations. This truncation treats quarks and gluons on an equal footing. 
It also takes into account the nonperturbative constraint that the total color charge be conserved and vanishing 
[22]. Moreover, the connection to the known Coulomb gauge heavy quark limit [23-25] can be established. The gap 
equations will be analyzed numerically for a particular nonperturbative input and various aspects will be explored. 
In the case of the gluon sector, it will be seen that a nontrivial dynamical mass scale emerges and which is connected 
to the nonperturbative renormalization of the theory. For the quark sector, both the dynamical chiral symmetry 
breaking and the heavy quark limit will be discussed. 

The paper is organized as follows. We begin in Sec. II by introducing the first order formalism, including the 
total color charge constraint that arises from the resolution of the temporal zero modes inherent to Coulomb gauge. 
The Dyson-Schwinger equations and the leading order truncation scheme will be presented in Sec. HI. The analytic 
development of the Dyson-Schwinger equations, in particular the reduction to the static equations and the heavy 
quark limit, is given in Sec. IV. Numerical results for both the gluonic and quark sectors appear in Sec. V. We close 
the paper with a summary and discussion. 
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II. FIRST ORDER FORMALISM 



To begin, let us consider the generating functional integral and action for QCD and review various aspects of the 
first order formalism in Coulomb gauge. The description initially follows closely in the spirit of Refs. [3, 15, 22, 26]. 
The generating functional is written as 



pJ,X,x\ = Jv^expl^tSgcD+t j dx\p%a% + J^-A% + x 



cQax "I" QaxXo 



(2.1) 



where T>^ generically represents the functional integral measure over all fields present, with sources p, J, x, X for tlie 
various gluon and quark fields (see below). In this section, only the source p will be relevant to the discussion and so 
we set J = X = X = 0- The (Minkowski space) QCD action reads 



Sqcd = J dxi^c 



i'J^Dqx + ij-Dx - m 



(2.2) 



In the above, {q) q^^ represents the (conjugate) quark field with fundamental color, spin and flavor indices collectively 
denoted with the index /3 and position argument denoted with subscript x. The Dirac 7-matrices obey the Clifford 
algebra {7^,7''} = '^g'^'^ with metric g^"^ = diag(l,— 1) (we explicitly extract all minus signs associated with the 
metric such that the components of spatial vectors such as x are written with subscripts, i.e., Xj). The temporal and 
spatial components of the covariant derivative in the fundamental color representation are 



Dox = dox - ig(Tx = dox - igty^T"-, 
Dx = Vx+ igAx = Vx + igA%T'' 



(2.3) 



where cr" and A"^ are the temporal and spatial components of the gluon field, respectively, and where the superscript 
a denotes the color index in the adjoint representation. The generators T" obey [r",T''] = if°''"^T'^, where the f"'^" 
are the structure constants and we use the normalization Tr[T"T^] = 5"^/2. The chromoelectric and chromomagnetic 
fields are written in terms of the gluon field as 



= -doxAt-Dfa^ 



1 



B» = V. X A« - X Al 

with the spatial component of the covariant derivative in the adjoint color representation given by 



D 



ab 



S'^'Vx - gr^'A] 



acb 



(2.4) 



(2.5) 



The QCD action is invariant under gauge transforms of the type 

a = UaU^ - ^{doU)U\ 

A ^ = UAU^ + (VC/) U\ 
q ^ q'^ = Uq 



(2.6) 



where Ux = exp{— iS^T"} is a spacctime dependent element of the SU{Nc) group. Because of this invariance, the 
functional integral (in the absence of sources) contains a divergence by virtue of the integration over the gauge group. 
Whilst this is in principle a global factor which can be absorbed into the normalization, when calculating Green's 
functions it leads to the necessity for fixing the gauge. This is typically achieved in the continuum formalism via the 
Faddeev-Popov technique which involves inserting the following identity into the functional integral: 



j Ve6 (f 1^]) Det [M'^^x, y)] , M'^^x, y) 



591 



F=0 



The expression F[<j, A\ determines the gauge condition {F = 0) and we choose Coulomb gauge: 

F[a,A] := W-A, 



(2.7) 



(2.8) 
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for which the Faddeev-Popov kernel reads 



M'^\x,y)r.-V-Df5{x-y). 



(2.9) 



There are caveats to the identity Eq. (2.7), namely that when the gauge fixing is incomplete, zero modes of the 
Faddecv-Popov operator will arise and one encounters the Gribov problem [1]. In Coulomb gauge there is a special 
case: temporal zero modes corresponding to time dependent but spatially independent gauge transforms [22], arising 
because the Faddeev-Popov operator involves only spatial diff'erential operators, and for which 



(2.10) 



Clearly, such zero eigenvalues (there are — 1 of them at each time xq) for the Faddeev-Popov operator automatically 
result in a vanishing functional determinant and hence invalidate the identity Eq. (2.7). In this work, we shall not 
consider the more general case of spatially dependent zero modes (which lead to the existence of Gribov copies in the 
usual sense). Following [22], we modify the original identity to 



1 = j 1)65 [f 



Det [W^ix, y) 



where VO explicitly excludes the temporal zero modes 9{xo) and 

[M-''{x,y)] = Det [M-»{x,y)] 



(2.11) 



(2.12) 



is the determinant with such zero modes removed. Replacing the source term (p) for the temporal gluon field in the 
functional integral (as mentioned, the other sources play no role in this section and are set to zero), the generating 
functional for Coulomb gauge QCD is thus written 



Z[p]= J (v-A^Det -V-D expl^iSgcD + i J dxp^al 



(2.13) 



To proceed, it is useful to convert to the first order formalism [3, 15, 22]. This is achieved by introducing an 
auxiliary vector field (tF), noting the following functional integral identity for the chromoelectric part of the action: 



dx 



(2.14) 



expjt j dx^El-El^ = j I?frexp|z j , 
The new field is then split into components using 

^ ' (2.15) 
changing variables tt ^ n — Vcj) and integrating out the Lagrange multiplier. The generating functional now reads 



const = j V(t>VTex.p^-i j (ia;r° (Vx-tt^ + V^C) j- 



Z[p]= j V^S S (v-7f) Det 



-V-D 



(2.16) 



with the action 



= J dxq„ 



S' = 



dx 



1 



1 



a0 



qi3x, 



2^x-B,-- 



^-.^i+^i.d^ai + \rx^lrx 



= j dxali^ 



]al3q0x^ ■ 



(2.17) 



Were the 5-functional constraints written in terms of Lagrange multiplicir fields, the determinant written in terms 
of ghosts and sources for all fields present, the theory would be in a local form suitable for discussing perturbation 
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theory as in Refs. [15, 18, 20]. The important feature about the decomposition to the first order formalism is that the 
action is now linear in a, which can be integrated out to give 



Z[p]= J VM(v -Ajslv-n^ Det -V-D S (v-D^ + -pj ei'^''+'^') 



(2.18) 



where 



(2.19) 



The 4> field can be integrated out by using the eigenfunctions of the Faddeev- Popov operator as a complete orthonormal 
basis for an expansion, the crucial point being that one must take into account the existence of the temporal zero 
modes. Following Ref. [22], the result is 



j V4>5 (y-D<k + p) exp |^ j dx \<kl^l4>l^ = S (^j dxp^ Det 



-V-D 



exp 



where 



(2.20) 



(2.21) 



The (5-functional constraint that emerges constrains the total color charge, the spatial integral arising from the 
projection onto the temporal zero mode. The generating functional is now. 



Z\p\ = j V^5 5 (v-7f) 5 (^J dxp 



with the action term 



/ 



dx 



2^x-^x +^x-"OxA^ - ^Px^x Px 



(2.22) 



(2.23) 



Some comments are in order. The (modified) Faddeev-Popov determinant in the generating functional cancels 
exactly, which is equivalent to saying that Coulomb gauge is ghost-frce and this arises from the elimination of the 
temporal gluon field. Notice though that the inverse Faddeev-Popov operator still plays a role. What is left of the gluon 
sector concerns the two spatially transverse vector fields A and if (tt would classically correspond to the momentum 
conjugate of A) which in quantum electrodynamics would give; rise to the two physical transverse polarization states 
of the photon. The remnant (5-functional constraint is the statement that the total color charge of the system must 
be conserved and vanishing and this explicitly includes a contribution from the external source p, the gluon field and 
the quark term. This is nothing other than the application of Gauss' law. 

In order to make sense of the generating functional, Eq. (2.22), the (5- functional constraint on the total charge must 
be rewritten in a useful form. This is most conveniently done with a Gaussian form: 



5(^1 dxpj ~^lim^AA(C)exp|-^ j dy dzr{y)C6''^S{yo - zq] 



W{z) 



(2.24) 



where C is a constant, M{C) is a normalization factor (that will be henceforth included implicitly in the functional 
integral measure) and the limit C — >■ 00 will be taken only at the end of any calculation. The generating functional 
can thus be written in the form 

Z[p]~^lim j P$5(v-i)<5(v-7r)exp05"}exp|-^ j dydzp^iv) (yo; ^) + C^"*"] %o - ^o)/(^)| (2.25) 



where 



S" =Sa+ dx 



■ 2 • -Dx - TT^^x -^x+T^x- OOxA^ 



=*a =*a I =;a 
— TT • TT ^ - 
2 



(2.26) 



and 



F"" {yo;y,z}^F^^'6{y-z). 



(2.27) 
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We shall refer to F as the Coulomb kernel. F is defined at eaeh time, jjo, since it involves no temporal operators. We 
see that the effect of the total charge conservation constraint is the presence of a potentially ill-defined and divergent 
constant additive term to the Coulomb kernel. It will be seen later that the constant, C, cancels out when considering 
physical quantities. 

As noted in Ref. [26], the Coulomb kernel is intimately related to the temporal gluon propagator. Although the 
temporal gluon field has been integrated out, the propagator still exists as a functional second derivative of Z[p] (this 
is why we have so far retained the source) and is defined in configuration space as 



red, ^_ 1 S^Z[p] 



Z[p] SipiSipd 



(2.28) 

p=0 



The temporal gluon propagator is known perturbatively in both the first and second order Coulomb gauge formalisms 
[18, 19]. Indeed, since the manipulations involved in going from the second to first order formalisms are integral 
identities performed on the generating functional, the spatial and temporal gluon propagators are identical in both 
formalisms. Taking two functional derivatives of Z given by Eq. (2.25) with respect to the source p and subsequently 
setting p = 0, one obtains: 



WfAv,w) = ^JvM (v-1) 5 (v-tt) exp{*5"} 



X hm exp <{ / dz p%y) [F"" (j/q; y, z) + Cc^""] 5{yo - z^)p\z) 

I {F''\v^; V, w) + CS"'^] 5{vo - wo) + J dy [F'^^vq; v,y}+ €6"'] p^vo, y) J dz [F'^^{wo; w,z)+ CS'^f] pf{wo,^^ 

(2.29) 

where 

p: = gr'^A'.-r, + gq^Al'nam^- (2-30) 

Were it not for the (5-functional constraint on the total charge, this expression would be the same as in Ref. [26] and the 
conclusion would be that the temporal gluon propagator splits into two parts, instantaneous and non-instantaneous: 

W^iiv,w) ^<iF'=''{vo;v,w)> S{vo - wo)+ <[F-p'=(t;o, «)] [F^^' pHwo,w)]>, (2.31) 

the instantaneous part being the expectation value of the Coulomb kernel and arising because the original Faddeev- 
Popov operator involves only spatial operators. However, the presence of the total charge constraint alters this: in 
particular, noting that the additional term is independent of the fields, we see that the instantaneous part of the 
temporal gluon propagator has the form 

W^iiv, wf^^ r.<iF'^{vo;v, w)> 5{vo - wq) + lim '£5'^5{vo - wq) (2.32) 

C— >cx) 

(in the original term involving the expectation value of F, the C — )• oo limit merely serves to reinstate the (5-functional 
charge constraint in the definition of the functional integral). The interpretation of the new term is simple: it is simply 
a (divergent) spatial constant and is completely nonperturbative in origin. As mentioned before, to make sense of the 
divergence, we shall consider C as being finite until the last step of any calculation. 



III. DYSON-SCHWINGER EQUATIONS AND TRUNCATION 

Let us discuss the Dyson-Schwinger equations and their truncation. The techniques involved in the derivation of 
the Dyson-Schwinger equations are standard, such that we shall present here only the most salient points in the 
interests of readability. The reader is referred to Refs. [15, 19, 20] for an explicit account of the derivation of such 
Dyson-Schwinger equations in Coulomb gauge. Noticing that the Coulomb kernel occurring in the action cannot be 
written as a finite order polynomial in the fields (due to the presence of the inverse Faddeev-Popov opcirator, this term 
is nonlocal), we must make some form of approximation in order to apply the standard Dyson-Schwinger formalism. 
To this end, we will introduce and justify a leading order truncation, whereby the Coulomb kernel occurring in the 
action is replaced by its expectation value and which will serve as an input into the resulting equations. 
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To start, let us consider the generating functional, Eq. (2.25), in the absence of sources. By implementing the 
transversality constraints on A, tt via Lagrange multiplier fields (A, r, respectively), the functional integral and corre- 
sponding action can be written 



-I' 



(3.1) 



where 



jx 



J dxdydz {-g [T'^ji]^^ (5(x - y)6{x - z)q^,A1yqp. + 5/"^^ [V ,J{x - z)] d{x - y)At,A%A^ 



2 fabc fade f Ab AC Ad Ae 
J J / ^ix^jx^ix^jx^ 

and where (recognizing the ^-dependence of the covariant derivative, D, given by Eq. (2.5)) 

F''''(^x,y;A) = [F'^' {xo; 5-,y) + CS^"] 6ixo - y,) 



(3.2) 



V X '-'x 



-Vx-D. 



cb 



6{x-y)+C6'''6ixo-yo). 



(3.3) 



Knowing that the Dyson- Schwinger equations are formed via functional derivatives of the generating functional, we 
observe that the A-dependence occurring within F always comes with an associated factor of the coupling, such that 
explicit functional derivatives of F will always result in additional loop structure. As a leading (loop) order Ansatz, 
we therefore make the following truncation for the Coulomb kernel: 



F'^" (x, y; a) ^ F'^" {x, y) = [F{x - y) + C] 5''H{x^ - yo) 



(3.4) 



where F is some scalar function which will serve as a nonperturbative input to the system of Dyson- Schwinger 
equations (the explicit expression will be discussed later). The Coulomb kernel is instantaneous (as is its expectation 
value) so that in momentum space, the function F will be independent of energy: 



F"^ {x,y)^J dke-'''<''-y^F''\k) ^ J d k e-'''<='-y^ 5'''' 



F{k^)+C{2^f5{k) 



(3.5) 



where dk ~ d"^ k / . It is clear from the discussion of the last section that F is intimately related to the 

instantaneous part of the dressed temporal gluon propagator. The tree-level contribution to F is ^ (~V^)~^ (or 

in momentum space), i.e., independent of A and so the one-loop perturbative results are in principle preserved within 

this Ansatz. 

Before continuing, it is worth briefly contrasting the formalism above with those of previous studies, namely 
Refs. [15, 20]. These studies focused on a local form of the Coulomb gauge first order formalism, whereby the a and ^ 
fields were not integrated out and the Faddccv-Popov determinant was written in terms of ghost fields. This form is 
ideal for studying perturbation theory. Here, after integrating out the a and fields, the Faddeev-Popov determinant 
cancels, leaving a form for the action involving the nonlocal Coulomb interaction term pFp. The purely spatial gluonic 
and quark components of the action (i.e., those that exclusively involve only the A, tt, q and q fields) are unaltered, 
meaning that many of the previous results from Refs. [15, 20] pertaining to these components are retained. After 
replacing F with its expectation value (which serves as an external input into the system), the Coulomb interaction 
term involves three new momentum dependent tree- level interactions between the A, tt and quark fields (their explicit 
forms will be presented shortly). These new interaction terms replace the dynamical interaction content of the a, 
(j) and ghost degrees of freedom with a simple set of effective vertices. In effect, the full nonperturbative towers of 
Dyson-Schwinger equations involving the cr, (j) and ghost fields have been 'solved' by the Ansatz for F. Obviously, this 
'solution' is only a leading order Ansatz; we shall however see that important nonperturbative physics is nonetheless 
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contained. As has been seen, the removal of the cr, (f> and ghost fields in the nonlocal formalism is intimately related 
to the imposal of Gauss' law and (total) charge conservation; in the local formulation, the Slavnov- Taylor identities 
perform this role [21]. 

The Dyson-Schwinger equations are integral equations relating the various proper (one-particle irreducible, [IPI]) 
Green's functions of the theory. The most basic of the IPI functions are the two-point proper functions and in the 
current formalism, with only A, n and quark fields (aside from the trivial Lagrange multiplier fields) , the most general 
momentum space decomposition of these is given by (see Refs. [15, 20] for details of the derivation and notation) 



^Tnrij (^) 



T^ab 

Airij 



(fc) = S'^^'ko 



Sij^V!r{k) + Zij(A;)r.n-7r(^) 

(k)+hj{k)TAAk) 
U,{k)TAA{k)+kj{k)TAA{k) 



<0) 
■ qqa0 



(k) 



j°koAt{k) - ^-kAsik) - Bm{k) + -f°koi-kAd{k) 



a/3 



(3.6) 



where lij{k) = kikj/J? is the longitudinal spatial projector and tij{k) = Sij — lij{k) is the transverse spatial projector. 
All (scalar and dimensionless with the exception of Bm and Ad) dressing functions are functions of k^ and }? separately 
due to the inherent noncovariance of Coulomb gauge. At tree-level, the dressing functions reduce to (the trcc-lcvcl 
forms for the proper two-point functions follow directly from the quadratic part of the action, 5^, given in Eq. (3.2)) 



A.. 



1, En 



Ad = 0. 



(3.7) 



Alongside the proper two-point hmctions, one is also interested in the corresponding propagators (connected two- 
point Green's functions). The connection between the connected and proper Green's functions is supplied via the 
Legendre transform. Since the components of the gluon field are treated individually in Coulomb gauge, the gluon 
propagator dressing function is not simply the inverse of the corresponding proper dressing function (as in Landau 
gauge) but rather, a matrix inversion structure arises. Taking into account the Lagrange multiplier fields that enforce 
the transversality of the A and 7f-fields, the components of the gluon propagator are given by (see also Ref. [15]) 



WfM^{k) 



iS^hijik) 



^inrik) 

'AM' 



Wf,,,{k) = t6'^''k%{k) 



A,(fc) 



where the common denominator factor, Ag(fc), including the Feynman prescription, is given by 

Ag(fc) = k^Tl^ik) - k^rAAik)T^^ik)+iO+. 
In the case of the quarks (see Ref. [20]) we have 



(3.8) 



(3.9) 



Wqqai3{k) = -- 



j\oAt{k)-j-kAs{k) + Bm{k)+j\o'y-kAa{k) 



A/(fc) 

where the denominator factor is given by 

Aj{k) = k^A'^ik) - k'Alik) - Bl{k) + klk'Alik) + ^0+. 



(3.10) 



(3.11) 



The matrix inversion structure for the gluon and quark propagators will have important consequences when solving the 
Dyson-Schwinger equations. At tree-level, the denominator structures reduce to the familiar forms: Ag = kQ — k^+iO+ 
and Af = kg — k"^ — + tO+. 

The tree-level vertex (three- and four-point proper) functions can be derived from the cubic and quartic parts 

of the action {S^ and 5^ of Eq. (3.2), respectively). As for the two-point functions, the purely spatial tree-level 
vertex functions of the local first order formalism can be taken from Refs. [15, 20]. For these, the momentum space 
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expressions are 

^AAMjkikuk2, ks) = -igr^" [5ij{k^ - k2)k + Sjk{k2 - k3)i + 5ki{ks - fci),] , 

^AAAMM^^^^^^^ ^3, ki) = w [SiMr^'f'"' - p'^-r^-) + dikSjiir'^f""' - p'^-f--) 

^^SL^i{kuk2, fcs) = -9 „^ , (3.12) 

where it is understood that all momenta are incoming and energy-momentum conservation has been applied (^ ki = 
0). The now interaction terms that arise in the present nonlocal formalism, with the Ansatz that the Coulomb kernel 
be replaced by its expectation value, are all contained within the quartic component of the action (<S** of Eq. (3.2)). 
They are all linear in F and since the color charge (p) involves both gluonic and quark components on the same 
footing, the trcc-lcvcl vertices all have the same structure. It will be seen later that indeed, the resulting Dyson- 
Schwinger equations for the gluon and quark sectors have very similar forms. Using the techniques of Refs. [15, 20], 
the momentum space expressions for the vertices are: 

^AA'l%kii^uk2, k3,ki) = -ig^ [r^ff'^dudjkF^fik, + k,) + r'^^'^'^jAkF^^iki + k^)] , 
4°«,(fci,^2,fc3,fc4) = ^g^r'^[^''Tf]^^5,,F'^f{k, + k2), 

^tLcMsi-kuk^MM) = -z5'[7°T«]^^[7°T'']^^F''*'(fci+M + *fl'[/T'']„,[7°r^],^^''"(A=i+fc4)- (3-13) 

Notice that the above vertices are written such that the symmetry properties are manifest. 

Let us now discuss the Dyson-Schwinger equations themselves. As stated earlier, we shall not present details 
of their derivation here: the basic techniques of such a derivation in Coulomb gauge are expounded in Refs. [15, 
20]. Gcncrically, the Dyson-Schwinger equations have a very definite structure of loop terms that arises from the 
combination of repeated functional differentiation of the generating functional and the Legendre transform connecting 
the connected and proper Green's functions. Once the notation and conventions have been established, different 
interaction terms merely follow this characteristic pattern; the difficulty is simply in keeping track of the various 
coefficients, especially where anticommuting Grassmann- valued fields (quarks or ghosts) are present. 

The Dyson-Schwinger equations for the two-point proper functions (in momentum space) are presented below. To 
aid presentation, the two-loop terms that will not explicitly be used in this study are omitted (they are collectively 
denoted F^^^ below). The equations are 

Tt,,{k) = Tf^^{k) - \ j d^T^^2t:l,,.{u:,-^,k,-k)Wi^^^^^^^ (3.14) 

+ j ^^4?i^a«i('^>-'^,-fc,fc)w%«/3«M + rt!i:^W. (3-15) 

^fAiAk) = ^fAiik)-l J ^u;ri°ltJ,,(fc,-fc,a;,-a;)W;f^,,(a;)-i J dujT^ZltMf''-k''''-'^)'^AAM 
+1 J ^AAMkiik, -oj, w - k)WXBkk- {^)Wrciv {k - ^)TUAVk'j{k - -k) 
- 1 ^^'^fX^ii^ - k, fc)H%,/3/3'('^)r|,A/3'a',-(^, k - uj, -k)W^,^,^{uj -k)+ T^^Mjik), (3.16) 

r^,a;3(fc) = r^^l^{k) + J dujr!^%^^^^,{k,-k,c.,-Lo)w-,,s,{uj) - 1 du;r^^X»m('''-'''^'-^)^'AM 

+ j ^^^fjLa'kik, -u>,w- k)Wgqa'0'HW^%kk'{k - u>)T^gB0,^k'{^, -k, k-u>) 

+r£a^(fc)- (3.17) 

Because of the existence of the mixed gluon propagator Wav and the possibility of dressed vertices such as F^^^ or 

TqqTt, in certain terms of the above equations (in particular, the one-loop term of Eq. (3.16) involving the three-gluon 
vertex and the spatial one-loop component of the quark self-energy in the above) one must sum up over the possible 



A 




FIG. 1: Dyson-Schwinger equations for F-^-^, F^a and Taa, omitting two-loop terms. Wavy lines denote proper functions, the 
large filled blob indicating the dressed function. Springs denote gluonic propagators, lines denote the quark propagator and all 
internal propagators are dressed. Small blobs indicate tree-level vertices and large circles denote dressed vertices. The gluonic 
field types B and C denote the sum over A and vf contributions arising due to the presence of mixed gluon propagators. See 
text for details. 



gluonic field types A and vf : this sum is denoted by the repeated subscript indices B,C , . . . (obviously, no relation to 
the chromomagnetic or ghost fields). The equations are presented diagrammatically in Figs. 1 and 2. As mentioned 
previously, the Ansatz of replacing the Coulomb kernel with its expectation value in the action does not interfere with 
the one-loop Dyson-Schwinger equations. Inserting the appropriate tree-level vertices and propagators into the above 
equations and using 

F°''(fc) = ,5'''' J- (3.18) 

(equivalent to the tree-level temporal gluon propagator) , it can indeed be verified that the known one- loop perturbative 
expressions [18, 20] are recovered. 

To complete this section, let us now describe the remainder of the truncation scheme. We are interested here in 
the leading loop, nonperturbative infrared behavior; in particular, on the effect of the terms generated directly by the 
Coulomb kernel, F. As discussed earlier, when replacing the full (nonlocal) Coulomb kernel with its expectation value 
within the action, one has already truncated out certain higher loop terms. The resulting Dyson-Schwinger equations 
retain their full leading (one-loop perturbative) structure. The next truncation we make on the Dyson-Schwinger 
equations is thus to restrict to the one-loop terms. 

The next part of the truncation scheme is to omit those terms generated by the tree- level r^]^^, r^]^^^ and rig^^ 
vertices. These are the terms that are present at one-loop in the Dyson-Schwinger equations, but do not involve 
the Coulomb kernel directly. What remains of the Dyson-Schwinger equations are the tadpole terms involving the 
tree-level four-point vertices T^j^^^, ri^^^^ and ^'^qqq'- these terms explicitly involve only the Coulomb kernel F and 
the propagators, thus forming a closed set of equations for a given input function F . 
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FIG. 2: Dyson-Schwinger equation for the quark two-point function, omitting two-loop terms. On the left-hand side, the filled 
blob indicates the dressed (inverse) propagator. Lines denote quark propagators, springs denote gluon propagators and all 
internal propagators are dressed. Small blobs indicate tree-level vertices and large circles denote dressed vertices. The gluonic 
field type B represents the sum over A and n contributions arising due to the presence of mixed gluon propagators. See text 
for details. 



In order to complete the truncation scheme for the Dyson-Schwinger equations, we must supply an expression for 
the Coulomb kernel. This is based on an infrared divergent 1/fc^ behavior. In this study, we shall not include the 
perturbative {1/k'^) term, since we are primarily interested in the infrared; however, an interesting point of note is 
that logarithmic ultraviolet divergences will still emerge (this will be discussed later). The 1/k* factor is justified on 
two grounds: lattice studies do show this behavior for the instantaneous component of the temporal gluon propagator, 
although the extant results admittedly do not go far enough into the infrared regime for this to be definitive [27-33]. 
The second reason is somewhat more pragmatic: as will be seen, such an infrared behavior is not only consistent with 
a linearly rising potential between heavy quarks but also generates dynamical chiral symmetry breaking for chiral 
quarks. Further, the product g'^F is a renormalization group invariant quantity [3, 26] and serves as an ideal quantity 
to use as input. For technical reasons, we use two forms for the input Coulomb kernel. The expressions are: 

g-^CpF^-^k) = (5"''C(27r)3,5(fc) +(5"''i^(fc2), 

f TTiTI (gluon/quark DSe), 

where the coefficients have been chosen for later convenience {Cp — {N^ — l)/2iVc is a color factor). The coefficient 
tJc is the Coulomb string tension (as distinct from the physical Wilson string tension, see e.g., Refs. [33, 34] for a 
discussion). In the Dyson-Schwinger equations, an unadulterated 1/k'^ factor results in an infrared divergence, so that 
we will introduce a numerical infrared cutoff regulator. In the case of the quark, a second infrared mass regulator, ^, 
will be considered — this form of regularization will turn out to have a very clear physical interpretation, which we 
will discuss later. Recall that the constant C multiplying the (5-function (arising from the total charge constraint) is 
considered to be finite until the end of the calculation whereupon we consider the limit C ^ oo. 



IV. TRUNCATED DYSON-SCHWINGER EQUATIONS: ANALYTIC DEVELOPMENT 



Let us now take the truncated Dyson-Schwinger equations and decompose them into a form useful for further 
analysis. 
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A. T-^A equation 
The truncated Dyson-Schwinger equation for F^a, Eq. (3.15), reads 

Inserting the tree- level vertex functions using Eq. (3.13) and the general forms for the propagators and proper two- 
point function, Eqs. (3.8,3.6), one obtains 



+g^r' [7°T^]„^ Sji J -^^fi^ - ^) [7°'^o^t(w) - l-ojA,{u:) + i?„(w) + 7°wo7-cDAd(w)]^„ . (4.2) 

We notice the occurrence of factors F[ijj — uj), which are technically undefined. However, the corresponding color 

structures read y^^^d^cd [T'^]q,^(5/3c, both of which are automatically zero. Assuming that F is infrared regulated, 
these terms would then pose no problem and we summarily dismiss them. Expanding the denominator factor Ag 
using Eq. (3.9), the equation then reads 

6,TA.ik) + uAk)VA.ik) = 6, - .,^iv. / ^(,) -XlHL(a.) + .0,] - ^'-'^ 

Since F is independent of the energy and the scalar dressing functions T aa, ^a-k^ ^Tt-n are all even functions of 
energy, the integrand is overall odd in the energy and thus vanishes without further assumption. We thus have the 
nonperturbative result that under the current truncation scheme, 

r^^(fc) = i, r^^(fc) = o. (4.4) 

In other words, the mixed gluon two-point proper function remains bare. This result is very useful because as will be 
seen below, the energy dependence of the gluon propagators turns out to be trivial. 

B. r,r7r and Vaa- the gluon gap equation 

The truncated Dyson-Schwinger equations for F^^ and Taa-, Eqs. (3.14,3.16), can be treated simultaneously. They 
read 

r«t,,(fc) = rW«^fc)-l| Au:Tf^^:l,,.{u^,-^,k,-k)Wi^Aik{^\ 

^tM,{k) = Tf;^{k)-\j du;rf^^:tjkiik,-k,u;,-co)W^^M^). (4.5) 

Notice that having restricted the truncation scheme to include only those one-loop terms involving the Coulomb 
kernel, the two equations are identical in their structure. Inserting the tree-level vertex functions, Eq. (3.13), the 
two-point function decompositions, Eqs. (3.6,3.8,3.9), projecting onto the transverse components (the longitudinal 
components of the proper two-point functions will play no role) and using the previous result that r^^r = 1 within 
this truncation scheme, one readily obtains 
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Given that F is energy independent, the energy integrals of the above are relatively trivial and since there is no 
/cQ-dependence, the proper dressing functions are energy independent. However, we should point out that really, the 
energy integral is only trivial if the spatial functions are regularized and finite. Here, this refers to an implicit infrared 
regularization in the case of _F ^ with its strong infrared singularity, the finite coefficient C multiplying S{k) and 
more generally, the ultraviolet cutoff when one considers the perturbative term F ~ 1/fc^ although we shall see that 
the Taa equation involves a logarithmic UV divergence even with the 1/k^ interaction. It is helpful to define the 
static propagators M^^'') as the energy integral of the full propagators (in configuration space, these are the equaltime 
propagators). In the case of the Waa propagator. 



witj{k) = zS"%jik) 



dkf) 



r^TT (^^) 



27r 



fc2-fc2^A^(fc2)^^,(/e2) + ^o+ 



= ^"%(fc) 



f 

J — ( 



dki 



27r 



kl + k''VAA{k^)T^^m 



(4.7) 

where in the second integral form, a Wick rotation (fco iki) has been performed. Doing the integral, one finds that 
the static propagator can be written in terms of a single dressing function, which we denote G: 



H41«;(fc-) = rt,,.(fc) 1 G(^^)V2, 



G{k^) 



Taa{P)' 



(4.8) 



The static propagator can also be written in terms of G: 



wit;{k) = s^''u,{ky4G{p)- 



1/2 



(4.9) 



The reduction of the two dressing functions Faa and T-^t^ to a single function G follows directly from the previous 
result Taw = 1- 

Returning to the Dyson-Schwinger equations, we can now write 



r,,(fc2) = l + ]g^N, 



d 



G{u'^Y'^F(k - uj)tji{k)Uj(uj), 



did 



I'U! 



2 fc2 



G(a;2)-^/2F(fc - uj)tji{k)tij{uj) 



(4.10) 



where duj = duj/{27r)^. Expanding out F with Eq. (3.19) we have 



r.,(fc2) = i+i^_^G(fc2)i/2 



+ 



TAAik') = l + l^-^G{Pr/' 



11% 

4Cf 



/ ^G{u^y^^F{k-co)tjiik)Uj{uj), 



\^ I l^^G{Q^)-''^F{k-a)t,S)UM- 



;2 A;2 



(4.11) 



The proper dressing functions, T-^jj and T aa (and subsequently, the corresponding propagators), arc explicitly depen- 
dent not only on C, but also include potentially infrared divergent contributions stemming from the integrals involving 
F. However, using the definition of G, one can easily see that the above equations can be written in terms of a single 
equation for G: 



G{P) = 1 



1 TV, 
4Gf 



- uj)tji{k)tij{uj) 



G{co^) 



2n1/2 



u;2 



G(fc2) 



fc2 G(d;2)i/2 



(4.12) 



The terms proportional to C cancel, showing that G and the static propagators are independent of C. In addition, 
one sees that the infrared divergence of F is tempered, such that G may be infrared finite (as will be seen). It 
would thus appear that the physical dynamics are contained within the static propagator dressing function, whereas 
the full propagators (and in particular, their pole positions) arc not physical. We shall discuss this at the end. For 
reasons that will become obvious shortly, we shall refer to Eq. (4.12) as the gluon gap equation. The form of the 
equation is identical to that for the static gluon propagator derived from the canonical approach [6] (see also the 
earlier work of Rcfs. [4, 13]). In the canonical approach, the Coulomb kernel F (here an input) is represented by a 
combination of ghost dressing and Coulomb form factors, which are self-consistently determined from their respective 
equations. Later work within the canonical approach included further terms: the ghost 'curvature' [7-9], and three- 
and foiu-gluon interactions [10]. The equivalence of the truncated Dyson-Schwinger equation above, Eq. (4.12), to its 
counterpart arising from the canonical formalism and considered in Ref. [6] is one of the results of this paper. 
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C. Tqq: the quark gap equation 

Under truncation, the Dyson-Schwinger equation for Tqq, Eq. (3.17), reads 

Inserting the appropriate tree-level vertices, Eq. (3.13) and propagators, Eqs. (3.8,3.10), resolving the color algebra 
(again discarding terms where one has Tr[r''] or f"''^ that multiply F{io — co) as previously discussed for Ta-k), one 
obtains 

r-q,Mk) = r£)„^(fc) + g'Cp J b°u;oMio) - 7°7-'^7%(w) +Bm+ T^^'cooAdH] , (4.14) 

where it is again recognized that F is independent of the energy. Projecting out the Dirac components in the 
decomposition for Egg, Eq. (3.6), one obtains four equations for the dressing functions: 



At{k) = l-ig Cf , . . X , 

J koAfiuj) 

du) k-ujAs(uj)F(k — w) 
k^^fico) ' 
Suj Bm{cj)F{k — Lj) 



A,{k) = l + tg^Cp 

Bm{k) = m + ig^Cp / 

Aa{k) = ig^Cp j 



2„ / duiu]Qk-LijAd(ijj)F{k — u)) 



kok^Af{cj) 



(4.15) 
(4.16) 
(4.17) 

(4.18) 



As for the Tatt equation, the odd energy integrals vanish since all dressing functions are functions of cOq (including 
A/), furnishing the result that 



At{k) = l, Ad{k)^0. 



(4.19) 



In addition, one sees that Ag and are independent of energy. The static quark propagator, defined in analogy to 
the gluon propagator, is then 



qqafi 



(fc) = -t[ 



dko 



j°ko-rkAgik^) + Bm{P 



al3 



2tt 



fc2-A;2A2(fc2)-B^(A;2) + zO+ 



7-M,(A:2)-B™(fc2) 



al3 



2n 



kl + k^Aj{k^)+Bl{k^) 



(4.20) 



and performing the integral, one obtains 



T,/tlt'2-\ _ ^rnik^) 



, M{k') 



2^P + M(fc2)2 ' Ag{P) 



(4.21) 



The static quark propagator can thus be written in terms of a single mass function, M. This is analogous to the case 
for the static gluon propagator. Expanding F with Eq. (3.19), the equations for Ag and Bm can be written as 



A,(fc2) = 1 



1 f dujk-LjF{k — Lj) 



I k- 



2 ^fc2 + M(fc2)2 2y fc2^^. + M(t«2)2 



Bm{k^) = m+- 



1 CM{P) 1 [ dujM(uj-)F(k -uj) 



+ - 



(4.22) 



Clearly the functions Ag, B^ (and hence the full quark propagator) are dependent on C and involve potentially 
infrared divergent integrals, just as for the gluon. Combining the above equations, we see that the C-dependence of 
M cancels: 



M{k^) = m + 



1 f dujF{k-Cj) 



2J y/uP + M{w^y 



(4.23) 
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and one again sees that although F may be strongly infrared divergent, M can still be infrared finite. Equation (4.23) 
is the quark gap equation and moreover, is well-known from the literature as the Adler-Davis truncation [12], where 
it was derived using the canonical Hamiltonian approach. A more advanced version of the quark gap equation to 
self-consistently include the spatial quark gluon vertex has been recently studied in the canonical approach [11]. The 
equivalence of Eq. (4.23) to the previously derived gap equation for the static propagator in the canonical approach 
[12] is again a result of this study. The similarity of the quark gap equation, Eq. (4.23), to its gluonic counterpart, 
Eq. (4.12), is striking (and this is why we take the liberty in referring to Eq. (4.12) as the gluon gap equation) and 
arises primarily from the truncation to include only those terms originating from the Coulomb kernel, which itself 
involves both the quark and gluon contributions to the color charge on equal footing. 

It is possible to make a connection between the (full) quark propagator and the leading order heavy quark propagator 
in Coulomb gauge [23] (see also [24, 25]). This is based on a spin-decomposition and we follow the spirit of Ref. [35]. 
Let us introduce the spin-projection operators 

P± = i(l±70), P++P_=l, P+P_=0, P|=P±. (4.24) 

These projectors furnish the following identities: 

P+7°P+ = P+P+, P_7°P_ = -P_P_, P+7°P_ = P+7*P+ = P_7*P_ = (4.25) 
which allow us to write the quark propagator as 

W-,,^p{k) = [{P+ + P_)W^,{k){P+ + P.)l^ 

(-^) 



x{\ko + Bm(fc^)l P+P+ - \ko - Sm(fc^)l P-P- - [P+7-fcP- + P-j-kP+\ As{k^)] . (4.26) 

LL JL JL J J a0 

The heavy quark limit for Coulomb gauge (implicitly in the rest frame) can be expressed as the limit |fc|/m — )■ 0. 
Considering Eq. (4.23), we can make a leading order estimate for the function M (and we shall see that this is 
confirmed by the numerical results). Given that F{k — 63) in the integral peaks at 63 = k, the infrared divergence is 
canceled and leaves 

M(fc^) « m + #^£l= l'^'" m + #. (4.27) 



+ M2(A:2) 

The functions Ag, given by Eq. (4.22) can also then be estimated: 

A(P) l + 0(l/m) 

B^(fc2) l'^^™ rn + ^C + ^J du3F{6j^) {= B^) (4.28) 

where in P^, it is recognized that the integral is infrared divergent and not suppressed by factors 1/m. We demand 
that \k\As — >• 0, despite the fact that C — >■ oo (in the heavy quark limit, m is the largest scale). The spin-decomposed 
quark propagator in the heavy quark limit is then 



[ko - Bra+ le] [fco + Brr,- Is] 

The first component represents a heavy quark propagating forward in time; the second, a heavy antiquark propagating 
backwards in time (and under time-reversal is equivalent to the first component). The first component is explicitly 
that found in Ref. [23] (and where only forward propagation is included) , showing that the leading loop order Dyson- 
Schwinger truncation scheme considered here reduces in the heavy quark limit to the scheme considered previously. 
Notice that when considering the heavy quark limit for the Bethe-Salpeter equation [23], the Faddeev equation [24] 
and the quark four-point Dyson- Schwinger equation [25], the constant C and the infrared divergence of the spatial 
integral over F that occur in B^, Eq. (4.28), cancel explicitly when considering color singlet quark configurations. 
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V. TRUNCATED DYSON-SCHWINGER EQUATIONS: NUMERICAL ANALYSIS 

Let us now consider the numerical solutions to the gap equations (4.12,4.23) for the static dressing functions G and 
M. The defining feature of the truncated gap equations is the input Coulomb kernel F 1/k'^. The strong singularity 
is also the defining problem in solving the equations. 



A. gluon gap equation 



In order to solve the gluon gap equation, Eq. (4.12), it proves convenient to change integration variables such that 
the radial integration momentum goes through the Coulomb kernel. The equation thus reads (to make the equation 
more readable, we use subscripts to denote the momentum dependence of the functions) 



Gk = l + 



4Cp 



--F^tji{k)tij{k - Q) 



G 



L/2 {k-Qf Gk 



k—u! 



P gT 

k—uj 



(5.1) 



'(fc-w)2 

The above equation is not renormalized nor is it regularized. It is convenient to use the following notation: 

2 



12 -'2 

k , y = Lo , 



{k-uj) = x + y- 2^/xyz, 



(47r)' 



dz 



(5.2) 



where both an ultraviolet (UV), A — >^ oo, and an infrared (IR), £ — >^ 0, spatial momentum cutoff have been introduced 
to regularize the integrals. Inserting the first form of Eq. (3.19) for F, the equation reads 



Gx = l + Ocl ] — 

V 



dy dz 



K(x,y;z) 



G. 



1/2 



Gx 



XQ 



1/2 



where 
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Kix y2) = l-fc^ 

7r(iV2-l)' ^^^'2^'^^ 26 • 



(5.3) 



(5.4) 



One immediately sees the effect of the strong IR singularity present in the kernel: as y — )■ 0, there must exist some 
cancellation such that the solution is independent of the IR-cutofi: (e). In effect, the combination of terms in the last 
bracket must cancel to leading order in y, such that the angular integral vanishes as y ^ (and for all values of x). 
Let us therefore consider the angular integral in some more detail: 



Iz{x,y) 



i: 



dz 



K{x,y;z) 



1/2 Gx 



For large values of x and vanishing y, 9 = x + 0{y) and Ge ^ Gx + 0{y). Then, it is obvious that 



Iz{x,y) 



dz 



1 ^/yx 



1 - 



2/(1-^') 



2x 



Gl'^0{y) 



(5.5) 



(5.6) 



which would then automatically lead to an overall IR convergent radial integral. The case where both x and y are 
small is slightly more complicated. In the process of regularization, an additional (mass) scale is introduced so let us 
assume that to leading order in the IR, 



where kq has dimension [mass]'^ and Go is some constant. The angular integral would then read 



Izix 



-^K{x,y;z) (-^) 



a/2 r 



i^y^'^oiy) 



(5.7) 



(5.8) 
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which can only have the necessary cancellation (given that x may be of the same order as y, but is still unfixed) for 
a = 1. Thus, the cancellation of the IR singularity leads us to consider a solution of the form 

= — ^ (5.9) 

where for low x, should be a constant. Indeed, in Ref. [6] where the gluon gap equation was solved with an infrared 

enhanced (but not l/k'^) Coulomb kernel, such a solution with a constant Kx was shown to be a good approximation 

1 /2 

to the full solution. Note that for such a solution, the static gluon propagator (~ Gx / ^/x) has a constant asymptotic 
value in the IR, characteristic of a massive solution and contradicting the analysis of [36] (see also, Ref. [37]), i.e., the 

propagator docs not have the Gribov form [1]. 

The gluon gap equation, Eq. (5.3) can be recast as an equation for Kx- 

Kx=c^c'y — dzK{x,y;zy _ , (5.10) 

Je y J-1 Vyv0 + Ks 

Let us discuss the UV behavior of Kx- Assuming that for finite, but large x, x ^ Kx, then 

~ crc7 / — / dzK{x, y; z) ^ . (5-11) 



y J-i Vy0 
For y ^ x (the upper limit of the radial integral), we would then have 



y J-1 2^ ' 3 
which after balancing dimensions leads us to the leading UV behavior for Kx- 

4 



ae7 / ^ ^ ^4(1 + z') ~ t'^cl ln(A) (5.12) 



4 / X \ 

~ -20-c7ln(^-^j . (5.13) 



Thus, although there is no perturbative content in the Coulomb kernel (~ ^/k^), there is still a logarithmic UV- 
divergent contribution to the static gluon propagator. Moreover, the factor 4/3, arising from the angular integral 
is characteristic to the ghost self-energy (see e.g., Ref. [38]). It appears rather ironic that the function Kx, which 
was introduced to cancel the leading IR singularity, is connected to the UV-cutoff, but in retrospect this is quite 
natural: the IR scale appears during regularization (breaking the scale invariance of the original theory). Further, the 
divergent term originates in the F^^-equation which, were one to have the perturbative form for F (^ 1/k'^), gives 
rise to a quadratic UV-divergence that is canceled by the gluon loop. Here, the quadratic divergence is reduced to a 
logarithmic divergence. One sees therefore that in Coulomb gauge, the nonperturbative UV-limit is not necessarily 
the same thing as the perturbative limit. 

Equation (5.10) is solved by iteration, to give the unrenormalized function Kx in terms of the UV-cutoff A. In order 
to numerically cope with the infrared divergence, the equation is rewritten in the form 

l + 72(a;;e,A) 

where 

dy [e-x + Kg] 



f dy f 

Ii{x;e,A) = ac-j — dzK{x,y;z) 

Je V J-1 

l2{x;e,A) = acj [ — I dzK{x,y;z) L (5.15) 

Je y J-1 y/yV^ + Kg 

The reason for using this form is that during iteration, where k within the integral is not yet the solution, the infrared 
cancellation required for Eq. (5.10) may not happen in practice and the error is amplified by the latent infrared 
singularity: the iteration procedure would thus be unstable. Both Ii and I2 diverge as e — > (the divergence is 
~ 1/ y/e in both cases) such that numerically their difference will be inherently prone to a loss of fidelity and the 
residues of the singularities are not known exactly enough to compensate. The ratio of the two integrals is better 
behaved such that with the form Eq. (5.14), the iteration is stable (but only slowly convergent). In the end, both 
forms of the equation are satisfied. 
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FIG. 3: Plot of for various values of the UV-cutoff, A. All dimensionful quantities are in units of ac- See text for details. 



A 


ao ao 


bi 


10*^ 


6.39 6.37 


-0.475 


10^ 


7.43 7.40 


-0.475 


10® 


8.47 8.44 


-0.474 


10" 


9.51 9.48 


-0.474 



TABLE I: Numerical values for fitted coefficients for various values of A. All dimensionful quantities are in units of ac- See 
text for details. 



As input, we take the values Nc = 3 and Cc = 1 (all dimensionful quantities are thus expressed in units of ac for 
now). The radial and angular integrals are performed using (standard) Gaussian quadrature grids. In the integrand, 
Kg is evaluated using either cubic spline interpolation derived from the grid of x-values or, for asymptotic arguments, 
extrapolation based on the following formulae: 



uq + aix + a2x'^ , small a; 
&o + ^1 ln(x), large x 



(5.16) 



In practice, the fit to the UV asymptotic form is performed for x € [A/10^,A/10] in order to avoid complications 
arising when x ~ A. It is found that for e < 10~^, becomes independent of e, as it should, is plotted in 
Fig. 3 for various values of A. It is seen that in the IR, goes to a constant value, whereas in the UV, it decreases 
logarithmically, confirming the previous analysis. The extracted values for the asymptotic coefficients ao (the infrared 
constant value of k^) and bi (the slope of the UV- logarithmic divergence) are given in Table I. The analytic value for 
bi according to Eq. (5.13) is —0.477 with the current input and one sees that this is numerically verified (to within a 
reasonable numerical precision of < 1%, recalling that the integrals themselves involve 0{l/y/e) contributions). Also 
given in Table I are values for an estimate for ao (ao)- This is generated by assuming that — ao is constant and 
then considering Eq. (5.10) for a; = 0: 



ao 



dy 



4 /4A 
o^cTln — 
3 \ ao 



This estimate is useful for understanding the A-dependence of the full solution Kx 
with A, but not quite; the same is true for k^. Empirically, we find that 

k(x') — k{x — x'aa{A); A) — ao(A) 



(5.17) 

ao scales almost logarithmically 

(5.18) 



is independent of A, shown in Fig. 4. 

The appearance of the scaled argument x' — x/ao{A) in Eq. (5.18) may at first seem somewhat arbitrary. We have 
seen that the presence of the ac/k^ term in the Coulomb kernel leads to the existence of an infrared mass scale ao 
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FIG. 4: Plot of k(x') for various values of the UV-cutoff, A. All dimensionful quantities are in units of cjc- See text for details. 



that is generated by the regularization scale A. The process of renormahzation is tantamount to choosing a reference 
scale, relative to which all other quantities are expressed. Importantly, both the function Kx and the string tension Uc 
have dimension [massY and must therefore be measured in the appropriate units. Choosing the (nonperturbatively 
generated) scale ao(A) as our reference scale (thereby implicitly choosing the renormahzation subtraction point /i = 
as will be seen below), we can rewrite the original equation for Eq. (5.10), in terms of rescaled quantities. Restoring 
the putative string tension, CTc, and denoting all scaled quantities with a prime: 

(5.19) 



ao(A)' " ao(A)' 
(similarly for y, 0, e and A itself), Eq. (5.10) can be rewritten as 



ao(A)' 



/«'(x';a^,A') = a^7 



A' 



V' J-i 



1 , 0'-x' + K'{e';<j',,A')-K'{x';a',,A') 
dzK[x , y ; z) 



(5.20) 



Notice that the equation for k! does not explicitly include reference to ao(A): the ao(A)-dependence resides in the 
condition 



ao(A) = k{x = 0; CTc, A) = ao(A)K'(a; = a;'ao(A) = 0; CTc = cr^ao(A), A = A'ao(A)) 



(5.21) 



(technically, oq is also dependent on CTc, but for notational convenience we shall drop the label) such that one may 
regard k' as a function of primed quantities and where k'(x' — 0; cr^. A') = 1. Subtracting at a;' = and defining 



AK'(a;', 0; o'^) = k'{x'- ct^. A') - «;'(0; a'^, A') = k'{x'- a^. A') - 1 
then we arrive at a closed expression for Ak'{x' , 0; u^): 



(5.22) 



l/y/1 + 0' + AK'{9\0;a'^) 



4 y' + AK'{y',0;a'^) 
3V7V^ + y' + A'^'(y',0;aO 



(5.23) 



Ak'{x' ,0]a'^) is explicitly independent of A', just as a renormalized quantity should be. It is however dependent on 
a'^, the scaled string tension. We see that, rather than in the unrenormalized case where everything is expressed in 
units of the string tension, in the renormalized case, everything is expressed in units of an implicit nonperturbatively 
generated scale. 



19 





2 









-2 


D 






-4 


"x 






-6 


< 






-8 




10 





=0, A'=10' 






=0, A'=10' 




— ii 


=0, A'=10" 




H 


=0, A'=l()" 






=(),A'=l(f, 


=0.5 


n 


=0,A'=10', a_. 


=0.25 


H 


=0, A'=10', a 


=0.1 


H 


= 10, A^IO' 






= 10, A^IO* 






10 



10"^ 10"^ 10° 10^ 



10 10 10 



FIG. 5: Ak'{x' , /i'; a'^) for various values of A', a'^ (= 1 unless otherwise stated) and two different fi' . See text for details. 



The subtraction at a; = a;' = is not the only possibility. Subtraction at x' = fi' is also possible and in this case 
the scale a(A) — k{x — /i'a(A); ctc, A) is chosen as a reference, such that 



and which obeys 



Ak'(x', At'; O = k'{x'; A') - ^'(/i'; A') = k'{x'; A') - 1 



A' 



A«'(x',/.';ai)=ai7 / ^ I J^zKix' ,y'; z) 



{x' ^ /.') 



(5.24) 



(5.25) 



The slope of the curves in the UV is proportional to a'^ and is the same for each value of /z' - the scaled string 
tension a'^ is independent of the renormalization point, although the scaling factor, a(A), is dependent on the choice 
of subtraction point /i'. This is nothing other than the statement that a'^ is the external input to the renormalized 
equation. Ak' {x' , fi' : a'^) is plotted in Fig. 5 for various values of A', a'^ and two different /i'. 

We notice that in Fig. 5, the curves for different fi' follow the same pattern as for the unrenormalized case. Fig. 3 
with different values of A. The following question arises: given Ak'{x' , /i'; cr^), can we derive the renormalized function 
for a different subtraction point, /I, and what are the associated scale factors? To answer the question, let us take a 
step back to the original unrenormalized function K{x;ac,A) and define 



a(A) = k{x = /i'a(A); ctc = o-^a(A), A = A'a(A)), 
6(A) = K(a; = 716(A); (Tc = CTc6(A), A = A6(A)). 



(5.26) 



So, knowing the unrenormalized function k{x; ac, a(A) and 6(A) can be derived for the chosen /x' and /i. Further 
defining the scaled functions as before: 



K(a:;ac,A) = a{A)K' {x — x' a{A); ac — <j'^a{A), A ~ A'a{A)), 
= b{A)K{x = xb{A);ac = CTc6(A), A = A6(A)), 

such that when written in terms of the appropriately scaled variables 

k'{x' — jj!; CT^, A') — 'K{x — JI; Wc, A) = 1. 

The renormalized (A-independent) functions are then 

AK'{x',fj,';<7'^) ^ k'{x ^ x'a{A);(Tc = (y'ca(A),A = A'a{A)) -I, 
A'K,(x,Ji;ac) = K(a; = x6(A);(Tc = CTc&(A),A = A6(A)) - 1. 



(5.27) 
(5.28) 

(5.29) 
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FIG. 6: Renormalized functions f{x) (labeled as jj,' = 0) and Ak{x, n;ac) (labeled /x = 6717) plotted versus x. See text for 
details. 

Importantly, both the above functions can be related to the original unrenormalized function k and the original 
variables through Eq. (5.27). One thus sees that 



Ak{x, /x; (Tc) + 1 



«(A) 
6(A) 



An' x' 



3(A) 
'a(A) 



(A); 



(5.30) 



So, one can indeed derive the renormalized function for subtraction point /x in terms of that subtracted at /x' and 
there is a single scale factor given by the ratio 



Z(^,/i';(7c, A) = 



a(A) _ k{x = ^'a(A); (Jc, A) 



(5.31) 



6(A) k{x = 716(A); (Jc, A) ' 

To test this, we take the unrenormalized function with CTc = 1, A = 10^ (plotted in Fig. 3) and consider the case 
/i' = to give a(A) = 8.4679 and cr^ = l/a(A) (implicitly in units of ctc). Considering /I = 6.7171 x 10'^, the ratio 
Z = a(A)/6(A) = 2.0067, with 6(A) = 4.2198 and = 1/6(A). The functions 



f{x) = Z [Ak\x' = x/Z, fi' = 0; cr^ = ajZ) + 1] - 1 

and Ak(x;/I;o'c) are plotted in Fig. 6 and one sees that they are identical (as they should be). 
Returning to the dressing function Gx, we see that 



(5.32) 



Gx = G{x;acA) 



x + K{x;ac,A) x' + 1 + An' {x' , ; a'^) 



G'{x',fi';a',). 



(5.33) 



When all dimensionful quantities are expressed in terms of the appropriate units (i.e., scaled by the nonperturbatively 
generated scale defined by the subtraction point fi'), the dressing function G is automatically independent of the 
UV-cutoff without modification and with a 'mass function' 1 + An' {x' , ^' ; cr'^) in these units. (The scaling factor 
Z(JI, 11'; ac, A), defined in terms of the unrenormalized solution k{x; CTcj A), tells us how the physical scale varies with 
different renormalization points.) In other words, the gluon dressing function requires no renormalization factor once 
the physical scale has been set. This has the implication that the nonperturbatively generated gluon mass would be an 
observable under the present truncation and with the 1 /q^ interaction, despite the naive appearance of a logarithmic 
UV-divergence. 
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B. quark gap equation 



As in the case for the gluon gap equation, it is convenient to change the momentum routing in the quark gap equation, 
Eq. (4.23), such that the radial integration momentum goes through the Coulomb kernel. The corresponding equation 
is then 



Mfe = m + 47r(T< 



/ 



1/2 



Mk- 



k ■ {k — uj) 



Mfe 



(5.34) 



where we have inserted the second form for the Coulomb kernel, Eq. (3.19): an infrared regularized form for the 

interaction with ^ having dimension [mass]^ and playing the role of a fictitious mass scale and for which the limit 
^ — >■ will be studied. (This is the infrared regularization method used, for example, in Ref. [39].) As has been 
emphasized, one can see the obvious similarities between the gluon and quark gap equations. In fact, there are some 
subtle differences arising from the specific forms for the integral kernel: in the quark case, the infrared singularities are 
somewhat more difficult to overcome than in the gluon case; however, the quark equation is explicitly UV-convergent 
with the above interaction. Notice that we will use the original form for the string tension, ctc, without mentioning 
the renormalization and scaling factors arising from the discussion of the gluon. For the quark case, all results may 
be expressed in units of ctc directly (as shall be seen) and one may use cTc or a'^ interchangeably as input. Using the 
conventions for the variables as for the gluon equation whereby P = x, etc., the equation reads 



M^=m + 



dz 



211/2 



Me- + J-zM:, 

X 



(5.35) 



In the above equation, a second possibility for infrared regularization emerges: one can set ^ = and study the 
infrared cutoff limit e ^ instead of the infrared mass regularization. It will be seen that both methods give identical 
results. 

The infrared analysis of Eq. (5.35) follows in the same manner as for the gluon. Setting ^ = for the moment it 
follows that for the radial integral to converge as y — )• 0, there must be some cancellation within the angular integral 
and we require that 



Iz{x,y) 



dz 



[e + Mi 



,1/2 



zM^ 



(5.36) 



vanishes fast enough in this limit, and for all values of x. Using the techniques as before, this is clearly achieved if 
Mx tends to a constant in the infrared (and as will be seen numerically). 

The UV analysis is best performed perturbatively. Consider the following integral (restoring here the original 
variables, fc, etc.): 



I = m + AwacMo 



k ■ uj 



W2]2 (fc-u;)2 + M2 



(5.37) 



which corresponds to the case of Eq. (5.34) where ^ = and where the function M within the integrand has been 
replaced by a constant, Mq (i.e., the above integral represents a first iteration of the full gap equation). The above 
integral can be performed using dimensional regularization (see, for example, Refs. [18, 20] for a discussion of such 
integrals) and the result for P ^ Mq is 



irk'^ 



(5.38) 



Importantly, the integral is explicitly UV-convergent (for all external momenta, k'^) and this means that the solution 
to the quark gap equation requires no renormalization, at least in the absence of the perturbative interaction ~ l/fc^. 
Unlike the gluon, the nonperturbatively generated dynamical scale thus plays no role in setting the physical units. In 
the full equation, Eq. (5.34), one would expect that the mass function vanishes in the UV like 1/k'^ when the bare 
quark mass m 7^ 0, since within the integrand, the mass function is negligible compared to the large momentum. In 
the chiral case, m = 0, one would expect that the mass function vanishes faster than l/P. 
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FIG. 7: Chiral quark mass function, M^, plotted for varying infrared cutoff regulator e. All dimensionful quantities are in 
appropriate units of the string tension, Gc- See text for details. 



Having discussed the asymptotic behavior, let us now turn to the numerical solution of Eq. (5.35). In the presence 
of the infrared singular integrals, it is again necessary to modify the equation for numerical use. To iterate, we use 
the form 



m + /i(a;;£,^, A) 



(5.39) 



where 



dzMg 



/2(x;e,$, A) 



[B + M\ 
dz 



211/2 



r2ll/2 



1 - 



(5.40) 



and we will consider two cases for the infrared regularization: e — !• with ^ = (infrared cutoff) and ^ — > with 
e ^ ^ (infrared mass). The iteration procedure for Eq. (5.39) is stable over a range of values for ^ and e, but becomes 
progressively more difficult as the regulators are made smaller and where the integrals diverge as l/y^ or 
The integrals are performed as for the gluon, the only difference being the asymptotic formula in the UV, where a 
powerlaw form 



(5.41) 



is used. Setting CTc = 1, all dimensionful quantities are henceforth expressed in units of CTc (or a'^ if one uses the 
renormalized scale defined via the massive gluon propagator). Note that the number of colors, iVc, has already been 
absorbed into the definition of the Coulomb kernel. Typically, the solutions to Eq. (5.39) are stable for A = 10'^ in 
the chiral case and A — 10^ for the massive case (the only exception is the heaviest mass m = 10, for which A = 10^ 
is necessary, see later for details). 

Starting with the chiral quark (to — 0) and with infrared cutoff regularization (i.e., studying the limit e — > with 
^ = 0), the quark mass function, M^, is plotted in Fig. 7 for various values of e. One sees that as e — 0, the mass 
function goes to a constant value M{x = 0) = Mq « 0.165 (in units of ^/oc) in the IR. The solution thus corresponds 
to a situation whereby chiral symmetry is dynamically broken. Numerically, the behavior in the UV corresponds to 
a powerlaw with exponent d ~ 1.84 (for the lowest value of e), which we shall discuss shortly. 

Repeating the analysis for the chiral quark, but with an infrared mass regularization (i.e., studying the limit ^ — > 0, 
with e < 10~® ^ ^), the solution is plotted in Fig. 8 for various values of ^. It is seen that as the regularization is 
removed, the mass function becomes identical to that using the infrared cutoff method, showing that the two methods 
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FIG. 8: Chiral quark mass function, Mx, plotted for varying infrared mass regulator ^. All dimensionful quantities are in 
appropriate units of the string tension, Gc- See text for details. 




FIG. 9: Chiral quark mass function at zero momentum, M{x — 0), plotted as a function of the infrared mass regulator ^. All 
dimensionful quantities are in appropriate units of the string tension, ac- See text for details. 



agree. The infrared constant value is the same as before Mq ss 0.165, although the UV exponent is a little larger: 
d « 1.90. Both sets of results for the chiral quark may be directly compared to those of Ref. [39] or, allowing for a 
factor Cf = 4/3 in the definition of CTc, to the results of Refs. [11, 12]. Turning to the UV exponent, Ref. [12] showed 
that one should expect d = 2 (in other words M ^ l/fc^ in the UV) . The numerical results here are somewhat lower 
(although within a reasonable numerical precision); however, given that the function is vanishing so rapidly, the effect 
of the UV tail is negligible. 

The behavior of the mass function for varying infrared mass regulator is rather interesting. In Fig. 9, the infrared 
constant mass, AIq, is plotted as a function of ^. It is seen that as ^ 0, Mq tends to its constant, nonzero value 
and indicating that chiral symmetry is dynamically broken in the presence of the strongly infrared enhanced Coulomb 
kernel interaction (and which one would naively expect). However, what is also seen is that for large Mq rapidly 
decreases and vanishes altogether for ^ > 0.02. As is well known, the interaction F ^ <Jc/k^ corresponds to a linearly 
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FIG. 10: Quark mass function, Mx, plotted for various bare quark masses. All dimensionful quantities are in appropriate units 
of the string tension, CTc. See text for details. 



rising potential with a coefficient given by Uc- Tlie explicit expression, appropriate to Coulomb gauge in the heavy 
quark limit and under the truncation to neglect pure Yang-Mills vertices reads [23] (temporarily reinstating the string 
tension Uc and all other constants): 



duo { C(27r)3j(tD) + 



if 



(1 



1 



(5.42) 



where r is a length scale. We shall shortly show that this form is indeed the correct expression for the truncation 
scheme considered here. The (5-function term proportional to C arises from the charge conservation term originating 
in resolving the temporal zero modes inherent to Coulomb gauge. This term does not contribute to the potential: 
such a term was in fact considered in Ref. [12], to the same effect. With the infrared mass regularization, the integral 
is 



V{t) 



1 - e 



> 
oo. 



(5.43) 



In the limit ^ — > and for finite length scale r, V{r), is a good approximation to the linearly rising potential, which 
gets better as ^ decreases. However, as ^ increases (whilst keeping r fixed), V{r) flattens to a constant and this 
constant decreases as ^ increases. The restoration of chiral symmetry for large ^ is now rather obvious - for large ^, 
the long-range part of the potential is no longer linearly-rising, but is constant. Thus, within the truncation scheme 
here, chiral symmetry is always broken by the existence of a pure linearly rising potential, the dynamically generated 
quark mass being measured in units of One cannot speak of a critical string tension in this respect (as compared 
to the concept of a critical coupling). However, the dynamical breaking of chiral symmetry is sensitive to the details 
of the potential when there is a flattening at large range, caused in this case by the large value of the infrared mass 
regulator. 

Let us finally discuss the case when the quarks have a nonzero bare mass. In the absence of UV-divergences, the 
bare mass requires no corrections due to renormalization. We again set CTc = 1 such that all dimensionful quantities 
are measured in the appropriate unit of the string tension. Using the infrared cutoff method, with e — 10^^ and 
A = 10* (except for the case m — 10, where A = 10^ is used), the solution to the gap equation for various quark 
masses is shown in Fig. 10 (the chiral quark is shown for comparison). One sees that in all cases, the quark mass 
function is constant in the infrared. In the UV, the exponent of the powerlaw lies in the range d = 1.02 — 1.07. This 
is comparable to the expected value d ~ I obtained earlier from the perturbative analysis, though as for the chiral 
quarks, the UV tail is not particularly important. In Fig. 11, we plot the deviation of the mass function from the 
bare quark mass: — m. Two features emerge from this plot. The first is that as m increases, the infrared constant 
value Mq — m extends further into the UV region (this is why A must be increased for the largest value of m, such 
that the solution converges properly). The second feature is that for increasing bare mass, m, the difference of the 
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FIG. 11: Difference of the quark mass function and the bare quark mass, Mx — rn, plotted for various bare quark masses. All 
dimensionful quantities are in appropriate units of the string tension, CTc. See text for details. 

infrared value Mq — m initially increases, but then turns over such that for heavy quarks, this difference becomes 
significantly smaller. In the heavy quark limit, one thus sees that — m is suppressed, at least in this truncation. 
Moreover, as the difference Mx — m vanishes when to — >■ oo, then the question of the ordering of the limits m — > oo and 
A — >■ cx) becomes irrelevant. One can thus see that the relation Eq. (4.27) is valid such that the connection between 
the truncated quark gap equation studied in this work and the heavy quark limit studied analytically in Refs. [23-25] 
(where contributions from pure Yang-Mills vertices were neglected) is established. This justifies the earlier use of 
Eq. (5.42). 

VI. SUMMARY, DISCUSSION AND CONCLUSIONS 

In this study, we have considered Coulomb gauge quantum chromodynamics within the first order formalism, 
under a leading order truncation and concentrating on the infrared behavior. In fixing to Coulomb gauge, particular 
attention must be paid to the temporal zero modes of the Faddeev- Popov operator. Within the first order formalism in 
Coulomb gauge, the temporal component of the gauge field can be exactly integrated out, leading to the cancellation 
of the Faddeev-Popov determinant; what remains of the gluon sector concerns two transverse field degrees of freedom 
A and tt . The resolution of the temporal zero modes further leads to the constraint that the total color charge be 
conserved and vanishing. This constraint was written in Gaussian form, resulting in a constant shift in the Coulomb 
kernel, proportional to C — >■ oo. 

Having integrated out the temporal component of the gluon field, the action is nonlocal due to the presence of the 
inverse Faddeev-Popov operator occurring in the Coulomb kernel. In order to derive the Dyson-Schwinger equations, 
a leading order truncation was introduced whereby the nonlocal Coulomb kernel was replaced by its expectation value 
in the form of an input Ansatz. This truncation led to the appearance of a set of new (momentum dependent) four- 
point interaction vertices. These vertices effectively replace the dynamical content of the tower of Dyson-Schwinger 
equations and Slavnov- Taylor identities involving the temporal, longitudinal and ghost degrees of freedom in the local 
formalism. The resulting Dyson-Schwinger equations were further truncated to include (and subsequently study) only 
those one-loop tadpole terms involving the nonperturbative part of the input Ansatz for the Coulomb kernel, which 
was taken to be strongly infrared enhanced (~ 1/k^). 

It was found that the truncated equation for the mixed gluonic proper two-point function, Tt^a, leads immediately 
to the result that this component is trivial. This in turn meant that in the subsequent gluonic Dyson-Schwinger 
equations for Ftttt and Taa, the energy dependence could be resolved. Moreover, the static gluon propagators could 
be easily identified in terms of a single dressing function, G. The proper two-point dressing functions were dependent 
on the constant C (arising from the charge constraint) and included potentially infrared divergent integral contributions 
(see below for a discussion). However, combining the two Dyson-Schwinger equations led to an equation for G alone: 
Eq. (4.12). This "gluon gap equation" for the static dressing function G was independent of C and (as was numerically 
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verified) had a finite solution, despite tlie strongly infrared singular interaction. An almost identical situation was seen 
for the quark sector: the static propagator was given in terms of a single dressing function, M , and a gap equation 
arose, Eq. (4.23), which was independent of C and whose solution was finite, despite the fact that the Dyson-Schwinger 
equations for the proper dressing functions were explicitly dependent on C and involved infrared divergent integrals. 
Importantly, the static gluon and quark gap equations were identical in form to their counterparts derived in the 
canonical Hamiltonian approach [6, 12]. This allows for an equivalence to be established between the respective 
truncation schemes and approximations employed within the two approaches. 

Using spin projectors, it was possible to consider the heavy quark limit and it was seen that the (full, not static) 
quark propagator reduces to the known result [23] within this truncation. This is important because in the heavy 
quark limit, the role of the infrared divergent integrals is understood: when studying physical, color singlet quantities, 
such infrared divergences cancel whereas for propagators, the pole position is shifted to infinity and reflecting the fact 
that infinite energy is required for an unphysical colored object to exist in isolation. This in tiirn explains the constant, 
C, stemming from the charge constraint: it is merely an additional, fully nonperturbative and constant contribution 
to be added to the infrared divergence. Given that the input Ansatz for the Coulomb kernel was shown to be directly 
related to the instantaneous part of the temporal gluon propagator, which in Ref. [23] was shown to be related to the 
quark-antiquark potential, C is simply a constant shift of this potential. In Ref. [12], it was demonstrated that an 
arbitrary, constant shift in the potential has no observable consequence (at least as far as the chiral quarks studied 
therein were concerned). Here, the charge constraint arising from the incompleteness of the gauge fixing (the temporal 
zero modes) results in exactly such an unobservable, constant shift (considering the limit C — )• oo) in the potential, 
leaving physical quantities untouched and naively prohibiting unphysical, color charged quantities from existing in 
isolation (in the sense that infinite energy is required to create such an object). 

For both the gluonic and quark sectors, it was numerically found that given the infrared enhanced input Ansatz, the 
static propagator dressing functions exhibited the presence of a nonperturbatively generated dynamical mass scale. 
Both functions were nontrivially finite and constant in the infrared. In the case of the gluon, despite the absence 
of the perturbative components, a logarithmic ultraviolet divergence emerged. However, expressing the gluon gap 
equation in terms of the nonperturbatively generated dynamical scale, the dressing function was seen to require no 
renormalization. This would suggest that the dynamical gluon mass would be observable within this leading order 
truncation insofar as the renormalization is concerned (i.e., if one were to ignore the previous discussion about the 
charge constraint and infrared divergences connected to the full propagator). The result is in disagreement with the 
Gribov-Zwanziger confinement scenario [1-3] and the results of Ref. [36] . This is presumably a shortcoming of the 
leading order truncation scheme utilized in this study. 

In the case of the static quark propagator in the chiral limit, the gap equation was solved using two infrared 
regularization procedures and their equivalence was numerically demonstrated. The results are identical to previous 
studies [11, 12, 39]. The nonperturbatively generated mass scale corresponds to dynamical chiral symmetry breaking, 
although one should point out that it is long known that quantitatively the resulting chiral condensate is too small [12] 
(this subject was tackled in Ref. [11]). Of interest however for the chiral quarks, was the infrared mass regularization, 
where the regularization parameter, ^, was large: in this case, the nonperturbatively generated scale decreased as ^ 
increased and eventually disappeared. Utilizing the connection between the input Ansatz for the interaction and the 
quark-antiquark potential, the observed restoration of chiral symmetry could be intuitively explained as a fiattening of 
the long-range part of the potential when the interaction is regularized. With heavy quarks, it was explicitly verified 
numerically that in this truncation, the heavy quark limit emerges naturally. 

Clearly, despite the fact that we retain only the leading order contributions and that the input Ansatz for the 
Coulomb kernel excludes the perturbative content of the theory, the Dyson-Schwinger equations of Coulomb gauge 
in the first order formalism represent a powerful tool to study nonperturbative quantum chromodynamics. Equally 
clearly, one can see that much is missing. For example, it is known from the canonical Hamilton approach that: the 
inclusion of the ghost loop ('curvature') in the deep infrared region [8, 9] and the triple-gluon vertex in the mid- 
momentum region [10] of the static gluon propagator have important roles; also that the spatial quark-gluon vertex 
is necessary to obtain a reasonable estimate for the chiral quark condensate [11]. Now, in the canonical Hamiltonian 
approach, the truncation scheme is defined by an initial Ansatz for the vacuum wavefunctional, thereafter it is a matter 
of computational effort to obtain results. In comparison, the Dyson-Schwinger equations can be derived completely 
(in the sense that, in principle, all loop terms can be written down) but must subsequently be truncated in order to 
furnish useful equations. Having made the connection between the two approaches at leading order here, it seems 
promising that beyond leading order such comparison may provide useful insights into both approaches. 

The heavy quark limit of the Dyson-Schwinger equations in Coulomb gauge also seems a very promising avenue for 
further study. From the lattice, it is known that the physical, Wilson string tension is not simply the coefficient of 
the infrared singularity in the instantaneous temporal ghion propagator (see, for example, Ref. [33]). The Coulomb 
string tension is analytically known to be larger than the Wilson string tension, encapsulated in the statement: "No 
confinement without Coulomb confinement." [34]. Combining the heavy quark limit of the Dyson-Schwinger equations 
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beyond the leading order truncation presented here with the Bethe-Salpeter equation, one may hope to study the 
quark- ant iquark potential quantitatively. Thereafter, phenomenological application to the hadron spectrum would be 
a realistic proposition. 
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